Role of two-flavor color superconductor pairing in a three-flavor Nambu— Jona-Lasinio 

model with axial anomaly 
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The phase diagram of strongly interacting matter is studied within a three-flavor Nambu-Jona- 
Lasinio model, which contains the coupling between chiral and diquark condensates through the 
axial anomaly. Our results show that it is essential to include the 2SC phase in the analysis. While 
this is expected for realistic strange quark masses, we find that even for equal up, down, and strange 
bare quark masses, 2SC pairing can be favored due to spontaneous flavor-symmetry breaking by the 
axial anomaly. This can lead to a rich phase structure, including BCS- and BEC-like 2SC and CFL 
phases and new endpoints. On the other hand, the low-temperature critical endpoint, which was 
found earlier in the same model without 2SC pairing, is almost removed from the phase diagram 
and cannot be reached from the low-density chirally broken phase without crossing a preceding 
first-order phase boundary. For physical quark masses no additional critical endpoint is found. 



I. INTRODUCTION 

The phase diagram of strongly interacting matter is 
studied with great effort, both experimentally and theo- 
retically. From direct observations we know that at low 
temperature and low chemical potential chiral symme- 
try is spontaneously broken and hadrons are the rele- 
vant degrees of freedom. There are strong indications 
from heavy-ion experiments that this is changed at high 
temperatures where quarks and gluons become the rel- 
evant degrees of freedom and a so-called quark-gluon 
plasma (QGP) is formed (see Ref. [l| for an overview). 
This is also confirmed by lattice simulations of quan- 
tum chromodynamics (QCD) at nonvanishing tempera- 
ture @, yj- These calculations indicate that, for physi- 
cal quark masses and vanishing chemical potential, the 
hadronic phase and the QGP are connected by a smooth 
crossover [J]. 

At high densities and low temperatures, on the other 
hand, strongly interacting matter is expected to be a 
color superconductor, where the quarks form Cooper 
pairs (for reviews, see [5N-111]) . This can rigorously 
be shown for asymptotically high densities by apply- 
ing weak-coupling techniques to QCD J12 - ll5j . Unfor- 
tunately, these methods fail at more "moderate" den- 
sities, which are of phcnomenological interest. More- 
over, because of the "sign problem" , the regime of low 
temperature and nonvanishing chemical potential can- 
not be studied within lattice QCD. However, in model 
calculations one typically finds that the hadronic phase 
is bordered by a first-order phase boundary in the low- 
temperature region [la . Il7j [ . The combination of this re- 
sult with the notion of the crossover at zero chemical 
potential then leads to the standard picture of the phase 
diagram where the first-order phase transition ends at a 
critical end point (CEP). The latter has attracted consid- 
erable attention, as it is potentially detectable in heavy- 
ion experiments |18j . 

It is possible, however, that the phase boundary of 
the hadronic phase contains more than one CEP. For in- 



stance, it was found within a Nambu- Jona-Lasinio (NJL) 
model that imposing charge neutrality or the inclusion 
of vector interactions weaken the first-order phase tran- 
sition and can lead to a second end point at low tem- 
peratures near the chemical potential axis |19l - l21j . Be- 
side these mechanisms, it has been shown in a Ginzburg- 
Landau (GL) analysis [22h24J that the interaction be- 
tween the chiral and diquark condensates, induced by the 
axial anomaly, can also lead to such a low-temperature 
endpoint. This result is rather general, but it strongly de- 
pends on the values of the GL coefficients, which could 
not be determined within this framework. It is there- 
fore very interesting that it has recently been confirmed 
explicitly within an NJL- model study [25(. Thereby the 
authors have included a C/,A(l)-synimetry breaking inter- 
action term which connects chiral and diquark conden- 
sates. This term plays a crucial role in the GL anal ysis 
but is usually neglected in the NJL model. In Ref. [25[ 
it was found that it indeed leads to the emergence of a 
second CEP if the coupling is sufficiently strong. 

However, the calculations have been performed under 
the simplifying assumption of three quark flavors with 
equal masses. In this case it appears natural that color 
superconducting quark matter is realized in the color- 
flavor locked (CFL) phase [26|, where up, down, and 
strange quarks are paired in a very symmetric way. The 
authors of Ref. J25J have therefore considered only one 
common diquark condensate and one common chiral con- 
densate in their analysis. The same is true for the GL 
analysis of Refs. |23 - l24| in the case of equal quark masses. 
Of course, the assumption of equal quark masses is quite 
unrealistic, leading to the question whether the second 
CEP can still be found in a scenario where the strange 
quark mass is significantly larger than the masses of the 
up and down quarks. In order to investigate this ques- 
tion, some generalizations of the model are necessary: 
When the quark masses are different, the chiral conden- 
sates will take different values for different flavors and 
also the diquark condensates will depend on the flavor 
content of the paired quarks. In particular, at large dif- 
ferences between strange and nonstrange quark masses, 



the CFL pairing will become unfavored [27 29] and a two- 
flavor color superconductor (2SC) [30, [3l| will emerge, 
where only up and down quarks are paired. 

In the following we will analyze the effect of these gen- 
eralizations on the phase structure. Thereby it was our 
original motivation to investigate the effect of realistic 
mass differences. However, to our surprise, we find that 
even for equal bare quark masses the 2SC phase can be 
favored due to the coupling of the chiral and diquark 
condensates by the axial anomaly. As a result we find 
the scenario of Ref. [25[ never to be favored, not even for 
equal quark masses. 

In the following, this will be discussed in detail. In 
Sec. Hi] we briefly introduce the model and the parameters 
before presenting the numerical results in Sec. IIIII and 
concluding in Sec. [IVJ 



II. THE MODEL 



A. The Lagrangian 



We adopt the N JL-type Lagrangian of Ref. [25( , 
C = q(ip -m + lm )q + tf + cf + tf® + C™ , (1) 



which is based on a frequently used Lagrangian (e.g., [32 
l35j|). extended by the interaction term £^J. It describes 
the dynamics of a quark field q with three color (r,g,b) 
and three flavor (u,d,s) degrees of freedom. The current 
quark masses enter through the diagonal mass matrix 
m = diagj (m u , rrid, m s ) and /i is the quark chemical po- 
tential. In the present analysis, we do not impose electric 
or color charge neutrality constraints. These constraints 
play an important role in compact stars and can lead 
to a rather complicated phase structure |33 - [3a |. On the 
other hand, they are less important in heavy-ion colli- 
sions, where the matter does not need to be locally neu- 
tral. As a first step, we therefore use a common chemical 
potential for all quarks, in order to keep the analysis sim- 
ple. 

The Lagrangian Eq. ((T|) includes a four-point interac- 
tion in the quark- ant iquark channel, 






a=0 



(2) 



and a four-point interaction in the quark-quark channel, 

3 

4 4) = H J2i (^UljCf) (q T Ci l5 til jq ) 

+ (qtiljCf) (q T CUljq) ] . (3) 

Here G and H are dimensionful coupling constants, 
C = i7 2 7° is the charge conjugation matrix, and r a 
are the Gell- Mann matrices in flavor space, extended by 
tq = ^/2/3IL/. The flavor and color structure of the 



quark-quark interaction is generated by the antisymmet- 
ric matrices 



h 



A 7 . 



ta = 

h = 



-n, 
-A 5 , 



1-3 = -^2 , 



(4) 



where A a are the Gell-Mann matrices in color space. 



(4) 



.(4) 



The four-point interaction terms C x and C d are sym- 
metric under U (3)^ x U (3) L transformations in flavor 
space. In order to break the U (1) axial symmetry we 
include the standard six-point interaction term [371 138( 

£f = -K {det/ [q(l + 75 )q]+det f [9(1-75)9]} , 

which can be related to instanton effects. This term con- 
nects three incoming right-handed fields with three out- 
going left-handed fields, and vice versa. In principle, it 
can therefore couple to three quark- antiquark channels 
(e.g., (uu)(dd)(ss)) as well as to one quark-antiquark 
channel together with one diquark and one anti-diquark 
(e.g., (ud)(ud)(ss)). However, in Hartree approximation, 
which will be employed below, C x only connects quark- 
antiquark condensates but no diquark condensates. Since 
on the other hand the coupling between quark-antiquark 
and diquark condensates was found to be important in 
the GL analysis of Refs. [Hi|, the authors of Ref. [H 
have introduced a second six-point interaction. In our 
notation, Eq. $§, this term can be written as 



K 



/.(e) _^_ v V 



i.j,k—l i 



{qUl k ^±l5)Cq 7 



q I C(l± l5 )t ] l k q)(q l (l± l5 )q J ) 



(6) 



where in the last factor the indices i and j refer to the fla- 
vor components. C ^ has the same symmetries as Eq. ([5]) 
but the quark fields are ordered in such a way that in 
Hartree approximation the diquark and quark-antiquark 
condensates are connected. This reordering of the quark 
fields can be realized via a Fierz transformation, 1 of the 
instanton vertex, relating C x to C I . In this way one 
can also relate the new coupling constant K' to K. How- 
ever, following Ref. [25(, we will treat K 1 as a free inde- 
pendent parameter. 



B. Mean- field approximation and thermodynamic 
potential 

We are working in mean-field approximation by intro- 
ducing the scalar diquark condensates 



Si = (q T Cj 5 tikq) 



(7) 



1 Useful references in this context arc Refs. [39H4 



and the scalar antiquark-quark condensates 
4>i = (QiQi)i i = 1,2,3. 



(8) 

(6) 



In Hartree approximation the six-point interaction C d 
Eq. ([S]), then is simplified to 



C 



(6) 

Xd MF 



K' 3 



-\si\ q % qi - s* q C'^Ukq^ 



+975 U k Cq T Si<j)i + 2 \si\ 2 4>i 



(9) 



One can identify the first term as a contribution to the 
effective quark masses, while the following two terms con- 
tribute to the anomalous selfenergy of the propagator in 
the color superconducting phase. The last term does not 
depend on the fields and gives a constant contribution. 

Adding the other terms from Eq. ([1]) and using Nambu- 
Gorkov bispinors ^ T = 1/V2 (q, Cq T ) the full mean- 
field Lagrangian can be written as 



C MF = 4 ,5-l$_V 

with the inverse dressed propagator 



(10) 



S~\p) 



i + M7o - M J2i=i A *75 uk i r[] 
-TLi A i^tih V-M70-M j ' 



Here M is the diagonal mass matrix of the constituent 
quark masses with the components 



K' 
Mi = nii - 4G<f>i + K\tijk\(t>j<t>k + -^-\ s i\ 2 

The off-diagonal elements of Eq. (jlip include 

( K ' 

Finally C includes a field- independent term 

3 3 



(12) 



(13) 



V = 2GJ2< 



AKfa 



E * 



i=l 



K' 
~2~' 



(14) 

With these ingredients the thermodynamic potential in 
mean-field approximation becomes 

Q(T, M ) = -T£ /^Xi T rln^ _1(lW "'^ 



(2ny 2 



T 



V. 



where the sum is over fermionic Matsubara frequencies. 
For equal diquark condensates, si = S2 — S3 = s, and 
equal chiral condensates, <p± — 4>2 — <j>3 = 4>, it agrees 
with the thermodynamic potential of Ref. [231, where 
only this limit was considered. 

The self-consistent mean-field solutions are given by 
the stationary points of Q with respect to the conden- 
sates, i.e., by the solutions of the gap equations 



dcf>i dsi 



0, 



1,2,3. 



(16) 



The phase diagram is then constructed by taking the 
solution with the highest pressure (i.e., the lowest value 
of Q,) at each point in the fJt-T plane. 



C. Parameters 

The model defined above contains eight parameters: 
the three bare quark masses (m u , m^, m s ), the two four- 
point coupling constants (G, H) , the two six-point cou- 
pling constants (K,K'), and a regularization parameter. 
For the latter we take a sharp three-momentum cut-off 
A. 

To fix these parameters, we basically follow the proce- 
dure of Ref. [25(. Starting point are the parameters of 
Ref. [H, 

m u = rrid — 5.5 MeV , 
m s = 140.7 MeV , 
A = 602.3 MeV , 
,2 



G= 1.835/A^, 
K = 12.36/A 5 , 



(17) 



which have been obtained by fitting meson masses and 
decay constants in vacuum. The value of the quark-quark 
coupling H, which cannot be determined from vacuum 
meson properties, is taken from Ref. [pj, 



H = l.l^/K 2 



(18) 



Finally, as already mentioned, we treat the coupling 
constant K' as a free parameter. We will often take 
K' = 4.2 K, because with this value a second endpoint 
was found in Ref. |25j |. 

The authors of Ref. [25| have only studied the case 
of equal quark masses. To that end, they reduced the 
bare strange quark mass to the value of the up and down 
quark masses and then re-adjusted the antiquark-quark 
coupling G, so that the dynamical up and down quark 
masses M u = Aid remain unchanged in vacuum (param- 
eter set II of Ref. (H): 



m.. 



m d = 5.5 MeV , G = 1.918/A 2 



(19) 



In this way the vacuum values of the chiral condensates 
4>u = 4>d and of the pion mass and decay constant remain 
unchanged as well. 

Since we are interested in the effect of the strange 
quark mass on the phase structure we generalize this pro- 
cedure to the case of arbitrary values of m s . Starting 
with equal quark masses, Eq. (|19[) . we will increase the 
bare strange quark mass up to the more realistic case, 
Eq. (|17p . always keeping the dynamical vacuum masses 
M u = Md fixed by adjusting the coupling constant G. 



III. RESULTS 

In this section we present our numerical results, which 
mainly consist of a series of phase diagrams with differ- 
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FIG. 1: The phase diagram in the fj,-T plane for K' — A.2K 
and equal bare quark masses m u — ma = m s — 5.5 MeV, 
only allowing for one common diquark condensate s and one 
common chiral condensate <j>. Thick (red) solid lines denote 
first order phase transitions, thin (green) solid lines second 
order phase transitions. The dotted (blue) line indicates the 
BEC-BCS crossover line defined by M(T, y) - tt, where M = 
M u =M d = M S . 



ent six-point couplings K' and different choices for the 
strange quark mass. We begin with the case of equal bare 
quark masses in Sec. 1111 Al before introducing larger val- 
ues of m s in Sec. IIIIBI 



A. Equal bare masses 

As a basis for our investigations, we consider equal 
quark masses, m u = rrid = m s = 5.5 MeV and allow for 
only one common chiral condensate <pi — 4>2 — </>3 and 
one common diquark si = 82 = S3- For K' — A.2K we 
then reproduce the phase diagram obtained in Ref. |25j . 
as shown in Fig. [TJ The normal phases, where all di- 
quark condensates vanish, are denoted by 'xSB' in re- 
gions with large chiral condensates and by 'NQ' in regions 
with small chiral condensates. At low chemical potential 
and high temperature these two regions are connected by 
a crossover, while at higher chemical potential and lower 
temperature they are separated by a first-order phase 
boundary. Without diquark condensates the latter would 
continue down to the chemical-potential axis. However, 
when diquark condensates are included, the CFL phase 
is favored in the lower right part of the phase diagram, 
replacing the NQ phase in that regime. The first-order 
chiral phase transition then does no longer go down to 
zero temperature but ends inside the CFL phase, 2 while 
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FIG. 2: The same as Fig. [T] but with independent diquark 
condensates si, S2, S3, and chiral condensates <f>\, 02, 4>3- 



the latter is bordered to the xSB phase by a second-order 
phase boundary. 

As discussed in Ref. [25], this particular phase struc- 
ture is a consequence of the interaction term C ^ . With- 
out this term, i.e., for K 1 = 0, the xSB-CFL phase tran- 
sition is first order and basically a continuation of the 
first-order %SB-NQ phase transition (see Fig. fU below). 
In particular, the chiral condensate drops considerably at 
the phase boundary and is therefore small in the entire 
CFL phase. This is different for sufficiently large values 
of K', as in Fig. [T] In this case the coupling between 
chiral and diquark condensates leads to the existence of 
strongly bound diquarks in the x^B phase, which even- 
tually condense and then form a Bose-Einstein conden- 
sate (BEC) in the lower-/i part of the CFL phase (see 
also Ref. [43| ). Hence, the xSB-CFL phase boundary is 
simply the condensation line of the diquarks and there- 
fore the transition is continuous. In particular, the chi- 
ral condensate is still large at the phase boundary and 
only decreases considerably at somewhat higher values 
of /i. In a small regime near the boundary to the NQ 
phase this happens discontinuously and is just the con- 
tinuation of the first-order phase transition between the 
xSB and the NQ phase. At most temperatures, how- 
ever, the decrease is continuous and closely related to 
a BEC-BCS crossover where the BEC-like CFL phase 
('CFLbec') gets converted into a BCS-like CFL phase 
('CFLbcs')- The corresponding crossover line is indi- 
cated in the figure as well. Following Ref. [25| we have 
defined it as the line where the in-medium constituent 
quark masses are equal to /1. Note that with this def- 
inition the crossover line does not exactly run into the 
endpoint of the first-order phase boundary. 



2 In this article we use the term 'chiral phase transition' whenever 
the transition is mainly characterized by a change of the chiral 



condensates (pi. Note, however, that in the CFL phase chiral 
symmetry is always broken by the diquark condensates. 



Next we allow the different diquark condensates Sj 
and the chiral condensates (j>i to vary independently for 
i = 1,2,3. This opens the possibility for the formation of 
a 2SC phase (si = s 2 = 0, s 3 ^ 0). Still keeping the bare 
quark masses equal, we obtain the phase diagram shown 
in Fig. [2j Whereas the normal-conducting phases NQ and 
xSB occupy approximately the same regions of the phase 
diagram as in Fig. [TJ we find that a large part of the CFL 
phase is replaced by a 2SC phase. The first-order chiral 
phase transition, which in Fig. Q] ends inside the CFL 
phase, now separates a BEC-like 2SC phase with large 
chiral condensates ('2SCbec') and a BCS-like 2SC phase 
with small chiral condensates ('2SCbcs')- However, un- 
like in Fig. [TJ the BEC-like and the BCS-like phases are 
always separated by a first-order phase transition and 
never connected by a crossover. 

Since Si and s 2 vanish in the 2SC phase whereas they 
are equal to S3 in the CFL phase, the two phases are nec- 
essarily separated by first-order phase boundaries. On 
the other hand, the phase transition from the CFL phase 
as well as from the neighboring region of the 2SC phase to 
the NQ phase is of second order. 3 As an interesting con- 
sequence, the first-order 2SC-CFL phase-transition line 
ends exactly on the second-order phase boundary to the 
NQ phase. Also note that the 2SC-NQ phase boundary, 
in the regime where it is second order, exactly agrees 
with the CFL-NQ phase boundary in Fig. [JJ This is due 
to the fact that on the second-order phase boundary all 
diquark condensates vanish and therefore the 2SC phase 
and the CFL phase have the same free energy. For similar 
reasons the second-order xSB-2SCbec phase boundary 
agrees with the xSB-CFLbec phase boundary in Fig. [1] 
Again, this is not surprising because both boundaries are 
related to the Bose-Einstein condensation of diquarks, 
which arc identical in the xSB phase and, hence, on the 
boundary. 

At low temperatures (T < 16 MeV) we find a small 
area between the 2SCbec phase and the 2SCbcs phase, 
where the CFLbec phase is slightly preferred. The origin 
of this "CFL island" will become more clear below. 

In order to get a more complete picture, we now vary 
the coupling constant K' of the six-point interaction £*, j . 
In Fig. [3] we show the phase diagram in the ji-K' plane for 
T = 0. One immediately recognizes that for all choices 
of K' there is always a first-order phase transition at 
some value of /1. For low values of K' only the xSB 
phase and the CFL phase are present. The first interest- 
ing development occurs when the interaction gets strong 
enough to have bound diquarks in the xSB phase, which 
then condense at some value of \x. As a consequence, 
a 2SCbec phase appears in the phase diagram. At a 
slightly higher value of K' the CFL phase splits into a 



^ 
k 




3 For simplicity, we drop the subscripts 'BCS' or 'BEC when the 
distinction between BCS-like or BEC-like phases is not relevant 
for the discussion. 



250 260 270 280 290 300 310 320 
u [MeV] 

FIG. 3: The phase diagram in the fi-K' plane at T — for 
equal quark masses, m u = rrid = rn 3 — 5.5 MeV. The mean- 
ing of the different line types is the same as in Fig. [TJ 



BEC-like and a BCS-like region. These two regimes are 
separated by a first-order transition, which eventually 
turns into a crossover at higher values of K' . The corre- 
sponding endpoint marks the critical value of K' where in 
the phase diagram with equal condensates (s± = s 2 = S3, 
4>i = §2 — ^3 ) the lower critical end point appears on 
the //-axis. However, allowing for 2SC pairing, there still 
is a first-order transition between the 2SCbec an d the 
CFLbec phase. 

Cutting Fig. E]at K' = 4.2 K, we reproduce the T = 
behavior of Fig. [5J From this perspective it becomes 
clear that the "CFL island" , which we have found there 
between the 2SCbec and the 2SCbcs phase, corresponds 
to the upper end of the CFLbec regime and will disap- 
pear at a slightly higher value of K' . 

In Fig. [4] we show a series of phase diagrams in the 
/i-T plane in order to illustrate how the phase structure 
evolves with K' . While the effect is almost negligible 
for small and moderate values of K', a major restructur- 
ing takes place at K' > 3 K: It starts with a 2SCbcs 
phase, which appears near the triple point of the three 
original phases (b) and then grows towards lower tem- 
peratures and higher chemical potentials. Slightly be- 
low K' = 3.5 if, Bose condensation of diquarks sets in, 
leading to a 2SCbec phase and shortly afterwards to a 
CFLbec phase (c). At the beginning, the latter is com- 
pletely separated from the CFLbcs phase by a first-order 
phase boundary, which turns into a crossover upon fur- 
ther increasing K' , As a result we find at K' = 3.744 K 
a critical endpoint (d) , which is of the same origin as the 
low-temperature critical endpoint in Fig. [T] Note, how- 
ever, that the endpoint is now located in a region which 
is separated from the low-density regime by a first-order 
phase transition between 2SC and CFL phase. More- 
over, the endpoint exists only in an extremely small in- 
terval of the coupling K': While for K' = 3.743 if the 
CFLbec-CFLbcs phase transition is still first order down 
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FIG. 4: The ^-T-phase diagram for different choices of K' and equal bare quark masses of (m u = rtid = rn 3 — 5.5 MeV). The 
meaning of the different line types is the same as in Fig. [l] In the 2SC phase the BEC-BCS crossover line (blue dotted line) is 
defined by the condition M Ut d{T, fi) = /i. 



to T — 0, already at K' — 3.745 K the entire phase 
boundary including the endpoint is covered by the 2SC 
phase. Eventually, the latter reaches the /z-axis and the 
two CFL regimes become separated again (e). At even 
higher K', the first-order 2SCbec-2SCbcs phase transi- 
tion turns into a crossover, leading to another endpoint 
(f). 4 This endpoint is located at the upper temperature 
end of the first-order boundary, like the "conventional" 
endpoint of the xSB-NQ phase boundary, and it moves 
downwards with increasing K'. Finally, slightly above 
K' = 5 K , it reaches the /^-axis and only a crossover re- 
mains (not shown). On the other hand, as argued before, 
the 2SC-CFL phase transition is necessarily first order, 
and therefore there is always a discontinuous phase tran- 
sition at low temperatures. 

After getting this overview, we would like to under- 
stand how the 2SC phase can be preferred over the CFL 
phase in some regions of the phase diagram although we 
have chosen equal bare quark masses. It clearly means 
that the Sf (3)-flavor symmetry is spontaneously bro- 
ken, so there is in fact a continuous set of degenerate 



4 As mentioned earlier, the BEC-BCS cross-over line (blue dotted 
line), which we have defined as the line where the non-strange 
constituent quark masses are equal to fi, does not exactly run 
into the endpoint. In this sense, the first-order phase transition 
(related to a discontinuity in the condensates) does not exactly 
coincide with the BEC-BCS transition, although both are closely 
related to each other. 



ground states, which are related to each other by flavor 
rotations. Bearing in mind that we will introduce larger 
strange quark masses later on, we choose, without loss of 
generality, the standard 2SC pairing pattern, i.e., S3 7^ 0, 
si = s 2 = 0. 

As obvious from Figs. [3] and 21 the emergence of the 
2SC phase at equal bare quark masses is a consequence 
of the six-point interaction C y ^. This term couples the 
chiral condensates to the diquark condensates and, thus, 
the constituent quark masses to the diquark gaps, see 
Eqs. (fT2"j) and (ITUl) . However, whereas in the CFL phase 
this happens symmetrically for all flavors, in the 2SC 
phase these equations conspire in a rather peculiar way: 
Since S3 7^ but si = S2 = 0, there is a contribution from 
the diquark condensates to M s but not to M u and M^. 
Hence, even if we start with equal bare quark masses, in 
the 2SC phase the strange quarks will be heavier than 
the non-strange quarks. In turn, the larger value of M s 
leads to an enhancement of A3 via an increased modulus 
of the (negative) condensate <j>3 in Eq. (|T^|) . 

This effect is illustrated in Fig. [5j where the diquark 
gap parameters and the constituent quark masses are 
shown as functions of the six-point coupling strength 
K 1 for the CFL and 2SC solutions at T = and \i = 
310 MeV. One can clearly see that A3 and M s in the 
2SC phase rise much faster with K' than the other quan- 
tities shown in the figure. In particular the ratio between 
A3 in the 2SC phase and the common gap parameter A 
in the CFL phase rises from about 1.25 at K' = to 2.0 
at K' = 4. It is therefore plausible that eventually the 
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stituent quark masses Mi (lower panel) in the 2SC and CFL 
solutions as functions of K' for m u — rrid = m s = 5.5 MeV 
at T — and /j, — 310 MeV. The vertical lines indicate the 
values of K' where the phase transitions from x^B to CFL or 
from CFL to 2SC take place (cf. Fig.©. 



gain in free energy is larger for 2SC pairing than for CFL 
pairing. In fact, if we had only diquark condensates and 
no dynamical quark masses, the pairing energy would 
be proportional to A| for each quasiparticle mode with 
pairing gap Aj. 2SC pairing would then be favored for 
A3I2SC > \/3A|cfl. Although in our case the dynamical 
quark masses (and their coupling to the diquark conden- 
sates) are crucial and certainly cannot be neglected, this 
estimate can at least roughly explain the numerical re- 
sults. 

At first sight the results found in this section seem 
to contradict the GL anlysis of Refs. 22-24] where no 
2SC-like solution was found for three equal flavors. How- 
ever, although starting from a general GL potential, the 
authors have chosen a particular ansatz for the conden- 
sates which does not allow for such solutions. By making 
a more general ansatz one can show that our results are 
indeed consistent with the GL analysis. This is discussed 
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in more detail in App. [X] 



B. Realistic strange quark mass 

We now introduce bare strange quark masses that are 
larger than the masses of the up and down quarks. In 
particular we are interested in m s = 140.7 MeV, which is 
the value obtained in Ref. |42| by fitting vacuum meson 
properties. We begin, however, with a more systematic 
investigation of the mass effects by gradually increasing 
m s from the equal-mass case m s = 5.5 MeV to the "re- 
alistic" value m s = 140.7 MeV. Thereby we adjust the 
coupling G for each value of m s as described in Sec. Ill CI 

The phase diagram in the p,-m s plane at T = for 
K' = 4.2 K is shown in Fig. |U As usual, we now de- 
fine the CFL phase as a phase where all three diquark 
gaps Aj take nonvanishing values but do not need to 
be equal. With rising m s , since pairs involving strange 
quarks become increasingly disfavored, the CFL phase 
gets more and more replaced by the 2SC phase, thus 
pushing the phase boundary to higher values of /1. As the 
spontaneous breaking of flavor SU(3) through the mech- 
anism described above is now stabilized by an explicit 
symmetry breaking, even a small enhancement of m s has 
a rather large effect. For similar reasons, the small CFL 
area between the 2SCbec phase and the 2SCbcs phase 
disappears already at m s w 7.5 MeV. Above this point 
there is a first-order phase transition between the two 
2SC phases, which ends at m s w 105 MeV. 

A larger value of m s also leads to a larger value of 
the strange chiral condensate </>3, which in turn leads to 
an enhancement of the effective four-point quark-quark 
coupling in the ud channel. The latter comes about from 
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closing a strange-quark loop in the six-point interaction 
term C ^ and adding this contribution to the genuine 

four-point vertex from C, . As a consequence, the di- 
quarks get bound more deeply, i.e., the diquark mass 
is reduced and, hence, their condensation, which deter- 
mines the transition from the %SB phase to the 2SCbec 
phase is shifted to lower quark chemical potential. 

In the remaining part of this article, we fix the strange 
quark mass at the "realistic" value, m s — 140.7 MeV. 
In Fig. [7] we show the phase structure at T = in de- 
pendence of the six-point coupling strength K' . With in- 
creasing K' the effective quark-quark coupling becomes 
stronger and the xSB-2SC phase transition moves to 
lower quark chemical potentials. Around K' = 2.2 K 
the 2SCbec phase forms at the low chemical potential 
side of the chiral phase transition line. With larger K' 
this region becomes broader and at K' — 4.0 K the first- 
order phase transition between the 2SCbec phase and 
the 2SCbcs phase ends. Similar to the equal- mass case 
we also find that an increase of K' is more favorable for 
the 2SC phase than for the CFL phase, so that the latter 
gets shifted to higher values of /i. 

In Fig. [HI we finally show a series of phase diagrams in 
the fi-T plane for different values of K'. A general result 
is that none of them contains a critical endpoint at the 
low-temperature side of the first-order chiral phase tran- 
sition. For K' = we recover the standard phase diagram 
(a), which has been calculated, e.g., in Ref. [29] for the 
same parameters: At low temperature chiral symmetry 
is restored in a first-order phase transition between xSB 
and 2SC phase, which is continued by a first-order phase 
transition between %SB and NQ phase at higher temper- 
atures and finally ends in a critical endpoint. Further- 
more, there is a CFL phase at low temperature and high 
chemical potential. With increasing K 1 the 2SC pairing 
becomes strengthened so that the CFL phase is pushed 



to higher chemical potentials while the first-order chiral 
phase transition gets successively "swallowed" by the ex- 
panding 2SC phase: As a first step the upwards-moving 
second-order 2SC-NQ phase boundary reaches the criti- 
cal endpoint of the xSB-NQ phase boundary (b). Next, a 
2SCbec phase emerges on the left-hand-side of the xSB- 
2SC phase boundary (c,d) and eventually the condensa- 
tion line unites with the 2SCbcs _ NQ phase boundary (e). 
At the same time first-order chiral symmetry restoration 
line becomes disconnected and ends inside the 2SC phase. 
When K' is further increased, the endpoint moves down- 
wards in temperature and finally disappears completely 



IV. CONCLUSIONS 

We have investigated the phase structure of strongly 
interacting matter within a three-flavor NJL-type model 
with an extended six-point interaction that couples chiral 
and diquark condensates. This interaction term, which is 
usually neglected in NJL-model studies of the phase dia- 
gram, has been introduced in Ref. [25| and mimics effects 
of the axial anomaly. While this in principle leads to a 
more complete picture, the corresponding coupling con- 
stant K' is basically unknown and was therefore treated 
as a free parameter in our analysis. 

Aim of the present study was to extend the investiga- 
tions of Ref. [25[ to realistic strange quark masses. In 
this context we have generalized the mean-field ansatz 
to allow for flavor-dependent chiral and diquark conden- 
sates. This opens the possibility for 2SC pairing, which 
was expected to become relevant at large strange quark 
masses. 

It turned out, however, that even for equal bare quark 
masses the 2SC phase is present in the phase diagram if 
the coupling K' is sufficiently strong. This spontaneous 
breaking of the SU(3) flavor symmetry occurs because 
the axial anomaly induces a mutual amplification of the 
strange chiral condensate and the non-strange diquark 
condensate. As an important consequence, we have not 
found a continuous transition from the low-density chi- 
rally broken phase to the CFL phase at T = for any 
value of K' between K' = and K' = 5 K. Related to 
this, the low-temperature critical end point which was 
found in Ref. [251 ] only survives in an extremely narrow 
K' interval and is otherwise covered by the 2SC phase. 

As expected, the 2SC phase plays an even more im- 
portant role at larger bare strange quark masses, in par- 
ticular for the "realistic value", taken from a fit to vac- 
uum observables. By varying K' we find several qualita- 
tively different phase diagrams, where the first-order chi- 
ral phase transition ends outside, inside or on the phase 
boundary of the 2SC phase or where there is no first-order 
chiral phase transition at all. It should be noted, how- 
ever, that most of these qualitative changes only occur at 
relatively large values of K' , which may turn out to be 
unrealistic. But even for smaller couplings, the anomaly 
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can be quantitatively important, as it stabilizes the 2SC 
pairing and shifts the CFL phase to higher chemical po- 
tentials. 



A. Acknowledgments 



The present study may be seen as a minimal extension 
of the analysis of Ref. 25] to include realistic strange 
quark masses. There are, however, many other aspects 
which have not yet been taken into account. Perhaps 
most important is the consideration of inhomogeneous 
phases. For a two-flavor NJL model without color super- 
conductivity it has been shown that the entire first-order 
chiral phase transition line is covered by an inhomoge- 
neous region and therefore removed from the phase di- 
agram [44h46I J . It would be interesting to see how this 
result is modified when strange quarks and color super- 
conducting phases are included and how the results de- 
pend on the axial anomaly. 



The work of H.B. was supported by the Hclmholtz In- 
ternational Center for FAIR and by the Helmholtz Grad- 
uate School for Hadron and Ion Research. M.B. acknowl- 
edges partial support by EMMI. 



Wc have also neglected the possibility of kaon conden- 
sation in the CFL phase 47j , which was found to be very 
important in an NJL model with realistic strange quark 
masses [3J, |35| . In these references, on the other hand, 

the anomaly effects described by C j have not been taken 
into account. This term should lead to a higher kaon 
mass and therefore suppress its condensation. It would 
be interesting to study this in more detail. 



Appendix A: Ginzburg-Landau approach 



The most general expression for the GL free energy up 
to order four is given in the equations (7), (13) and (14) 
in Ref. [23J as the difference to the normal phase 



Finally, we should recall that we have restricted our- 
selves to a single quark chemical potential. Similar in- 
vestigations of anomaly effects should also be performed 
for electrically and color neutral matter. 



n ($, d L , d R ) = n x ($) + n d (d R , d L ) + n xd ($, d R , d L ) 

(Al) 
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with 



il x =^Tr[$t$] + h -l (Tr[$t$]) 2 + ^Tr[($t$) 2 



4! 



Co 



(det[$] + det[$ t ]) , 



fid =a Tr[d L d' ! L + d R d R ] 

+ (3 1 ((Tr{d L dl}) 2 + (Tr{d R dl)y 



7t \2l 



+ f3 2 Ti[(d L dl) 2 + (d R d R ) 

+ p 3 Tr[(d R dl)(d L d R )] + l3 4 Tr[d L d[] Tr[d R d 

n xd =7iTr[(d i? 4)$ + (d L d R )&], 

+ A!Tr[(d L d[)$$t + ( dfl 4)$t$] j 

+ A 2 Tr[di4 + d R d R ] • Tr[$t$] 

+ A 3 (det[$] • Tr[(d L d R )$- 1 } + h.c. 



R\ ' 



(A2) 



Here $, dh and c£r are 3x3 matrices containing the 
chiral and the left- and right-handed diquark fields, re- 
spectively (for details, see Ref. [23J). This ansatz has 13 
unkown coefficients which are in general T and /i depen- 
dent functions and cannot be determined within the GL 
analysis. 

For the investigation of three massless flavors the au- 
thors of Ref. [23[ have therefore restricted themselves to 
an ansatz with equal condensates for all flavors, $ = 
diag (cr, <7, a) and di, = —d R = diag (d, d, d) . The free 
energy then takes the simplified form 



Cl 3F (a, d) 



—a H — cr 
3 4 



d 2 + ^d 4 ) -jd 2 a + \d 2 a 2 



(A3) 



with only seven independent coefficients, 5 which are re- 
lated to the original GL coefficients in Eq. (JA2I) by 

a = 3a , c = 3c , b = - (3&i + b 2 ) , 

a = 12a , = 12 (60i + 20 2 + 3 + 30 4 ) , 

7 = 671, A = 6(A 1 +3A 2 -A 3 ) . (A4) 

Yet, the full exploration of the remaining seven- 
dimensional parameter space would still be almost im- 
possible. However, by making additional, physically mo- 
tivated, assumptions (e.g., positivity of c, 0, 7, and 
A) interesting non-trivial results have been extracted in 
Refs. [HH. 

Similarly, two-flavor systems have been investigated 
in Ref. [23| choosing $ = diag (cr, cr, 0) and dh = —d R = 



diag (0,0, d), i.e., by completely neglecting all conden- 
sates which involve strange quarks. The resulting free 
energy is given by 



9l 2F (cr, d) 



r 7 + i ff 



^ + ^) +X 'd 2 a 2 , 
2 4 ' 



(A5) 



with five independent coefficients, which are different 
from those in Eq. (|A3j) . In this context it is impor- 
tant to note that by neglecting the strange condensates 
completely, it was assumed that the strange quarks are 
infinitely heavy and therefore decouple from the non- 
strange sector. In particular, a direct comparison of Q 2 p 
with Q3F in order to study the competion between 2SC 
and CFL pairing was not intended and is not possible. 

In contrast, in our NJL-model anlysis, we have stud- 
ied a 2SC phase in a three-flavor environment with a 
large strange quark chiral condensate. The correspond- 
ing GL ansatz is $ = diag (cr, a, a s ) and d^ — —d R = 
diag (0, 0, d), yielding 



9. 



2SC 

a 2cr 2 

2 I 



{a,a s ,d) 



—a cr, 
3 



K 



2 ^4 
— cr 



&i - 26 2 



^22 
cr a . 



f4* 



±d 2 cr s + X[d 2 a 2 + X' 2 d 2 a 2 , (A6) 



where the coefficients a, c, a, and 7 are the same as in 
Eq. (|A"4|) . and 



b' 2 = \{b 1 +b 2 ) 



6i = -(26 1 + 6 3 ), 

0' = 4 (2A + 20 2 + 3 + 4 ) 

Al = 2 (2A 2 - A 3 ) , A 2 = 2 (A! + A 2 ) 



(A7) 

Hence, instead of five GL coefficients as in Eq. (|A5|i . 
^2SC depends on nine independent coefficients. In par- 
ticular, unlike fl 2 p, it contains the anomaly induced c- 
and 7-terms which are related to the six-point interac- 
tions Eqs. ([5]) and (HJ) in the NJL-Lagrangian. 

Obviously, the extra terms in Eq. (|A6|) are due to the 
strange chiral condensate cr s , which is not present in 9t 2 p- 
We can thus recover Eq. (|A5|) by setting a s equal to zero. 
In fact, a more appropriate way to describe the decou- 
pling of the strange quarks with large masses is to treat 
a s as a nonvanishing constant. The 7-term, among oth- 
ers, then leads to a renormalization of the ^-coefficient 
a, while, e.g., the 6^-term leads to a shift of the vacuum 
energy. This underlines again that a diract comparison 
of the free energies obtained with Q 2 p and Q 3 p would be 
meaningless. 

On the other hand, keeping a s as a dynamical variable, 
a comparison of Q3F with £l 2 sc is possible. 6 To that end, 



5 In addition, a term of order <r 6 must be introduced by hand in 
order to stabilize the system if b < 0. For simplicity, we neglect 
such terms in the present qualitative discussion. 



6 Here we tacitly assume that the GL anaysis is meaningful in both 
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FIG. 9: The Ginzburg-Landau free energies of CFL solutions 
(solid) and 2SC solutions (dashed) as functions of 7 at the 
endpoint found in Ref. [23 - [24| for the example described in 
the text. 
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gies in the examples given in Fig. [9] 



7, X[, A 2 ), plus possible six-order terms which might be 
needed for stability reasons. An exhaustive GL analysis 
is therefore practically impossible. 

However, since we only want to show that 2SC-like 
solutions can indeed be favored in some cases, it is suf- 
ficient to demonstrate this by an example. In order to 
proceed, we therefore make a few simplifying assump- 
tions: (i) We neglect all A-terms, i.e, Ai = A2 = A = 0. 
This was also done in Ref. [23| in the three-flavor case, 
(ii) In the NJL model without the anomaly terms and 
without diquark coupling, i.e., taking into account only 
the interaction term C x , the different flavors do not mix. 
Hence their contributions to the thermodynamic poten- 
tial are additive. To reproduce this feature we choose 
b[ = 26/3, b 2 = 6/3. (iii) In the diquark sector we choose 
j3' = /3/(3-2 2//3 ). In the absence of chiral condensates (or 
for 7 = 0) this leads to the relation djsc — 2 1 / 3 g?3f, in 
agreement with weak-coupling QCD Q. As can be seen 
in Fig. O this relation is also approximately fulfilled in 
the NJL model at K' = 0. (iv) We restrict the possible 
2SC solutions to BCS-likc solutions, i.e., we set a = in 
Eq. ([Mil . 

For the 2SC phase, we then obtain the simplified GL 
potential 



^2SC {a s ,d) 

a 9 6 

-at H c 

6 s 12 



6 12 • 2 2 / 3 



^^2 

-a ere 



(A8) 



while the CFL phase is described by Eq. (|A3[) with A = 0. 
Thereby both potentials depend on the same coefficients 
so that the resulting free energies can be compared with 
each other. In order to further reduce the number of pa- 
rameters, we perform this comparison at the location of 
the low-temperature critical end point found in Ref. 122 
0, 



c 

36 



2f_ 



a 



/3c 2 
' 27 7 6 2 



(A9) 



we should replace the coefficient 6 in Eq. (|A3j) by the 
combination 3(6^ — b' 2 ) and the coefficient A by 3(A^ + 
A2). Moreover, we have to take into account that the 
coefficients /3 and /?' are linear independent (see Eqs. (|A4j) 
and (|A7|) ). This can be accounted for by writing j3 = 
3(3/31 + (3' 2 ) and (3' = (3[ + /3' 2 with p[ = 4(2/3i + /3 4 ) 
and /?2 = 4(2/3 2 + Pz)- Thus, in oder to compare the free 
energies of the 2SC and CFL solutions, we have to deal 
with ten independent parameters 7 (a, c, b[, 6' 2 , a, /3[, f3' 2 , 



phases, i.e., that all condensates are small. Clearly, this does not 
need to be the case. 
7 A straightforward way to obtain a unified description of both, 
CFL and 2SC phase, is to make the ansatz $ = diag (a, a, cr a ) 
and (Ll = —dn = diag (d! , d! , d) . In this case there arc eleven 



This eliminates a and a from the equations, so that we 
are left with four parameters, 6, c, /3, and 7. 

Finally, we introduce an arbitrary scaling factor A of 
dimension energy and measure all dimensionful quanti- 
ties in units of the corresponding power of A. We then 
simply choose 6 = c = /3 = lin these units and study 
the different free-energy solutions at the point given by 
Eq. (|A9|) for varying values of 7. This anaysis has been 
done numerically. 

In Fig. |H]we show the results for &2SC and Q3F for our 
example. We see that the CFL solution is favored at low 
values of 7 whereas the 2SC solution becomes favored at 



independent coefficients. Our approach is more restrictive and 
therefore has one parameter less. 
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larger values of 7. This phase transition is in qualitative 
agreement with our findings in the NJL model when K' 
is increased, cf. Fig. [3l 

The values of the condensates associated with in min- 
ima of the free energy are shown in Fig. 1101 Their be- 
havior with rising 7 is qualitatively similar to the effect 
of an increased K' in the NJL calculations, cf. Fig. [3] In 



both cases when increasing the coupling (7 or K') the 
diquark and the strange quark chiral condensate in the 
2SC phase are rising much faster than the condensates 
in the CFL phase. 

The NJL-model results are thus completely consistent 
with a GL analysis if the latter is performed with a suf- 
ficiently general ansatz for the condensates. 
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